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We introduce a one-dimensional (ID) XZ model with alternating 
of of +1 and ofof +1 interactions on even/odd bonds, interpolating between 
the Ising model and the quantum compass model. We present two ways of 
its exact solution by: (i) mapping to the quantum Ising models, and (ii) us- 
ing fermions with spin 1/2. In certain cases the nearest neighbor pseudospin 
correlations change discontinuously at the quantum phase transition, where 
one finds highly degenerate ground state of the ID compass model. 

PACS numbers: 75.10.Jm, 03.67.-a, 05.70.Fh, 63.70.Tg 

Introduction — Models of magnetism with exotic interactions are motivated 
by rather complex orbital superexchange in Mott insulators. In certain case the 
degeneracy of 3d orbitals is only partly lifted and the remaining orbital degrees of 
freedom are frequently described as 1/2 spins. They also arise from spin-orbital 
superexchange, with rich dynamics leading to enhanced quantum fluctuations 
near quantum phase transitions [1], and to entangled spin-orbital ground states 
[2]. The orbital interactions have much lower symmetry than the SU(2) of spin 
interactions and their form depends on the orientation of the bond in real space 
[3], so they may lead to orbital liquid in three dimensions [4]. A generic and 
simplest model of this type is so-called compass model introduced in [5] when the 
coupling along a given bond is Ising-like, but different spin components are active 
along particular bonds, for instance J x o-faj and J 2 af cj along a and b axis in the 
two-dimensional (2D) compass model [6]. This situation is quite different from 
classical Ising-like models with periodically distributed frustration [7]. 

Recently we discussed the properties of the one-dimensional (ID) model 
that interpolates between Ising and compass model [8]. The solution based on 
specific choice of interactions indicated divergences in correlation functions while 
approaching the transition point. This suggested first order phase transition but it 
was not clear if the effect is generic or it is just an artefact of singular parametriza- 
tion of interactions. To answer this question we introduce more general solution 
which can be applied to any parametrization and then we make second insight 
into the original problem. 
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XZ model in one dimension — The model is described by a generalized XZ 
Hamiltonian with different energies for even and odd bonds, where the number of 
sites is 2N and for simplicity we suppose that N is even 

N 

K = X] { Jl °2i-l cr 2i + J202i-l°2i + L ^2i a 2i+1 + L ^ u 2i u 2i+l\ ■ ( l ) 
i=\ 

This Hamiltonian turns into the one discussed in Ref. [8] if we fix the energy 
constants as follows: J\(a) = |(|1 — a\ + 1 — a), J2(a) = —J\(a) + 1, L\(a) = 
-Ji(-a + 2) + 1, L 2 (a) = Ji(-a + 2), where < a < 2. For < a < 1 it gives 

N 

= S {(1 - <x)02i-l°- Z 2i + + «+l} - (2) 

i=l 

and for 1 < a < 2 we transform erf of , {2i — 1, 2i} <-> {2i, 2i + 1} for all i and 
a -> (2 - a). 

firsf solution — To solve the model given by Eq. (2) we choose eigenbasis 
of of operators consisting of vectors \si, s 2 , S3, ... , S2iv) with Sj = ±1 for all z. 
Every state like that can be denoted equivalently as 

\ti,t2,... 1 t N ) rir2 rN = \t 1 ,t 1 r ll t 2l t 2 r2,...,t N ,t N r N ) , (3) 

where U = S2i-i and Ti = s 2 i-iS2i for i = 1, 2, . . . , N. This will let us exploit 
the fact that the hamiltonian (2) flips only odd pairs of spins. For states like (3), 
we define new spin operators T t z and rf which act only on t\, t2, ■ ■ ■ , ijv quantum 
numbers 

T f 1*1 j *2, • ■ • ! tff) rir2 ... rN — I — t\, t2, ■ ■ ■ , tjy) rir2 ... rN , 

T l Z \tl, *2, ■ ■ ■ , tjv)rir 2 --T N = *1 1*15*2, ••• ,tN)nr 2 ---r N ■ (4) 

Next we transform each r- z as follows r- 2 = rir 2 • • • rj_irf , to get effective forms 
of the Hamiltonian in subspaces spanned by vectors (3) with fixed 7*j 's 

AT 

H rir2 ...r N {a) = {«+! + a^f} + , (5) 

1=1 

where Tjy +1 = rir 2 . . . r^rf and Cp(a) = (1 — a) X^i=i r i- Now we have to solve 
two types of quantum Ising model (QIM); either with periodic or antipcriodic 
boundary condition. The solution is well known and was described in detail in 
Ref. [8]. First step is to introduce Jordan- Wigner (JW) transformation, 

Tj = (c J -+ct)JJ(l-2ct Ci ) ) 

i<j 

rf = (l-2c} Cj ). (6) 

The boundary condition for fermion operators {ci} differs for even and odd num- 
ber of quasiparticles in the chain. Fortunately, the Hamiltonian conserves the 
parity of fermions, this lets us write Hp = |(1 + P)T~i? + |(1 — P)T~L? , where 

JV 

ft$ = Ei( c l - c «)( c l+i - Ci+i) - 2acjci> + C ?( a ) > ( 7 ) 

i=l 
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and where |(1±P) are projections on subspaces with even (+) and odd (— ) num- 
ber of quasiparticles. Here we adopt notation r for the subspace labels rir 2 . . . rjy. 
The boundary conditions are cn+i = T c i YliLi r i f° r T~($ respectively. After the 
Fourier transformation, Cj — -^^2^^^°^ wc nn d (?) in a block diagonal form 

Hf = Y, r± {clc k {2 cos k - 2a) + (c{cl k e lk + h.c.)} + C?(a) + Na , (8) 
k 

where quasimomenta k take "integer" values k — 0, ±7^, ±2^, . . . ,7r for a pe- 
riodic boundary condition and in the antiperiodic case they are "half-integer" 
k = 0, ±-^, ±2-^, . . . ,7r. Diagonalization is completed by a Bogoliubov transfor- 
mation 7^ = MfcC^ + VfcC_fc preformed for all fc > and k ^ tt where (uk,Vk) are 
eigenmodes of the Bogoliubov-de Gennes equation [Wi,7t] = Ek"ft- In this way 
we get the full Hamiltonian's spectrum in every subspace f. For instance, in case 
when n»=i ri = 1 we obtain a Hamiltonian for even number of quasiparticles 7^ 

«r («) = Y, r+E * (^7* - + ( 9 ) 

where _E fc = 2{1 + a 2 — 2acosfc} 1 / 2 is the quasiparticle energy. The ground 
state is Bogoliubov vacuum in a subspace where all = — 1, apart from points 
a = 0, 1, 2 it has no degeneracy. Cases of a = 1,2 are trivial. For a ~ 1 we have 
Cf(a) — for all r which means that there is a ground state in every subspace 
where YiiLi r « = 1- This results in 2 7V ~ 1 -fold degeneracy for the ID compass 
model [Fig. 1(a)]. In the limit N — > 00 the lowest energies of periodic and 
antiperiodic QIM get equal at a = 1. For < a < 1 they are already two- and 
threefold degenerate, so when a = 1 the total degeneracy is 5 x 2 N , and the spin 
gap vanishes [8]. 

Second solution — The most direct way of dealing with Eq. (1) is to leave 
the interactions and Li, 2 undefined and to start with the JW transformation 

°j = (^ + ^(1-24^), 

i<j 

°*j = ](c j -c])Y[(l-2c\c i ), (10) 

i<j 

which transforms spins into fermion operators Cj. Here the crucial step is to 
introduce new quantum number for fermions with two possible values n and p. 
This can be regarded as quasiparticles' spin or as splitting the chain into bi- atomic 
elementary cells. We define c™ = cn-\ and c\ = c-n. Because of the boundary 
conditions and the fact that TL preserves the parity of fermions, we obtain 

N 

n± ^J2 { J+ -M + L + <?M + i + J-$W + L-cfMU + h-c) , (11) 
1=1 

where 7i + {Ti~) is defined in subspace with even (odd) number of fermions with 
c n+i = ~ c i ( c n+i = c i)- H ere we introduced new notation for J's; J* = 
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(a) (b) 

Fig. 1. (a) Eigenenergies E n of the XX-ZZ model (2) for 2N = 8 sites and for increasing 
a. Level crossing at a — 1 marks the quantum critical point of the ID compass model, 
(b) Intersite pseudospin correlations on odd {2i — 1, 2i} and even {2i, 2i + l} bonds in the 
XX-ZZ model (2) for increasing a. A transition between two types of pseudo-order, 
with (o'2i-iO'2i) = —1 for a < 1 and (p^i^i^-Ll — ~ 1 f° r ce > 1, occurs at the quantum 
critical point a = 1, where only {(J^i-i^i) — (ffjO'li+i) = — § are finite. These results 
are reproduced from Ref . [8] . 



n,p 



Ji ± J2 and analogically for L's. Now, we proceed with Fourier transform c-' 
^^X)fe ±e *"' fec fe' P l° r n an d P fermions separately, compatible to fermions' parity. 
For periodic boundary conditions k's take "integer" values and in the antiperiodic 
case they are "half-integer" . Finally we find the problem block diagonal in a form 



TL 



k 



+ {c n k^l(,r 



L- e - lk 



)) + c^4(J++L+e-> k ) + h.c] , (12) 



and for TL similarly but with k — 0, 7r possible. Diagonalization is completed by 
a four-dimensional Bogoliubov transformation. We search for invariant subspace 
of a linear operator [TL + , .] in 8-dimensional space spanned by c±^) and • The 
result suggests the form of transformation as 



( ^ \ ( C * \ 



»t 



7_fe 
pt 
% 



\7! 



= Pk 



< 1, 



(0 < k < tt), 



y»t 

ITT 
v Pt 



= Ar 



,(13) 



\CL k J 



where [3k and (3^ are orthogonal matrices. This assures that 7fc's are fermionic. 
The rows of j3 are eigenvectors of Bogoliubov-de Gennes equation [H,J^ ] = 
^klb and eigenvalues are energies of sytem's elementary excitations. Finally, 
we find TL + in a diagonal form 



pt p 



(14) 



where £™ = 2{ Jf 



Li 



2JiL 2 cosfc} 1 / 2 , = 2{jf + Lf + 2J 2 Licosfc} 1 / 2 . In 



a similar way we get the result for H . Luckily, the parity of particles is the 
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same as parity of original JW fermions, thus only states with even (odd) number 
of quasiparticles 7^ belong to the spectrum of TL + (H~). The ground state energy 
E obtained from (14) is E = -^J2 k + i E k + E l\- Putting J 1>2 = Ji, 2 («) and 
Li t 2 = Li^((x) we find the same energy spectrum as described in Ref. [8] [see 
Fig. 1(a)].' For a < 1 one finds E% = 2(1 - a) and E p k = {1 + a 2 - 2acosfc} 1 / 2 , 
which means that the occupation numbers 7^ 7^ play the role of the subspace 
indexes Ti from the previous solution while 7^7fc describe excitations within a 
given subspace. 

Pseudospin correlation functions can be derived from E as derivatives with 
respect to Ji, 2 and £1,2, respectively. The main result is that (cr| i _ 1 cr| i ) and 
{ a 2i a 2i+\) remain constant in intervals [0,1) and (1,2], but with discontinuities 
at a = 1 [see Fig. 1(b)]. The origin of these singularities is a cusp of E^{J\,L 2 ) 
surface at (J\,L 2 ) — (0,0). The trajectory (Ji,L 2 )(a) passes through this point 
at a = 1. This means that no first order phase transition occurs in the general 
model of Eq. (1), unless curves (Ji,L 2 )(a) or (J 2 ,L\){a) pass through (0,0). 
Passing should be interpreted literally as passing, not reversing at (0,0). For 
example, the curve (a;, y) = ((a — l) 2 , (a — l) 3 ) reaches the point (0, 0) at a = 1, 
but nevertheless one of the pseudospin correlations remains continuous. 

Summary — We have presented an exact solution which demonstrated that 
a hidden order with constant pseudospin correlations exists in the ID XZ model. 
The second method used for solving this problem provides more insight into the 
nature of the quantum phase transition as discussed in Ref. [8] , while the second 
one is more flexible and may be generalized, for instance, to the ladder geometry. 
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